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ABSTRACT 

Two  complementary  methods  of  describing  the  high-speed  rotor  noise 
problem  are  discussed.  The  first  method  uses  the  second-order  transonic 
potential  equation  to  define  and  characterize  the  nature  of  the  aerodynamic 
and  acoustic  fields  and  to  explain  the  appearance  of  radiating  shock  waves. 
The  second  employs  the  Ffowcs  Williams  and  Hawkings  equation  to  success¬ 
fully  calculate  the  acoustic  far-fleld.  Good  agreement  between  theoretical 
and  experimental  waveforms  Is  shown  for  transonic  hover  tip  Mach  numbers 
from  0.8  to  0.9. 


INTRODUCTK 


It  Is  now  well  known  that  linear  acoustic  theory  falls  to  correlate 
with  many  features  of  measured  transonic  rotor  noise.  In  particular,  linear 
theory  does  not  predict  measured  acoustic  pressure  levels  or  waveforms  of 
rotating  blades  operating  with  transonic  flow  In  the  tip  region  (Refs.  1-5). 
The  discrepancy,  In  level  and  waveform,  is  particularly  large  for  helicopter 
rotors  that  use  relatively  thick  airfoil  sections  (as  compared  with  pro¬ 
pellers)  operating  at  high  Mach  numbers.  As  a  result,  large  transonic  dis¬ 
turbance  fields  alter  and  enhance  the  radiated  acoustic  signature,  making 
high-speed  rotor  compressibility  noise  an  important  and  practical  research 
problem. 

Until  recent  years,  the  theoretical  approaches  of  rotor  noise  uti¬ 
lized  the  "acoustic  analogy"  of  Lighthill  (Ref.  6).  Ffowcs  Williams  and 
Hawkings  (Ref.  7)  extended  the  analogy  to  allow  for  relative  motion  between 
the  acoustic  sources  and  an  observer,  yielding  an  integral  equation  that 
describes  the  resulting  noise.  This  equation  was  refined  for  eventual 
numerical  computation  by  many  authors,  most  notably  those  of  References  8 
and  9.  Unfortunately,  the  helicopter  rotor  problem  Is  so  complex  geometri¬ 
cally  that  an  exact  "acoustic  analogy"  was  not  Initially  attempted.  Instead, 
the  governing  equations  were  linearized  and  the  resulting  Integrals  were 
numerically  evaluated. 

Although  the  success  of  the  linear  model  at  low  Mach  numbers  can 
still  be  debated.  Its  success  at  transonic  tip  Mach  numbers  is  quite 
limited  (see  Fig.  1) .  Neither  the  amplitude  nor  the  shape  of  the  far- 
fleld  acoustic  disturbance  is  accurately  predicted.  In  fact,  although  a 
radiating  discontinuity  (shock)  was  experimentally  measured,  linear  theory 
cannot,  to  any  extent,  account  for  this  phenomenon.  Similar  discrepancies 
have  been  confirmed  by  tests  with  both  model  and  full-scale  helicopters  In 
high-speed  forward  flight  (Refs.  1,  10). 


22-1 


Figure  1.  Comparison  of  theory  and  experi¬ 
mental  pressure-time  history,  In-plane, 
r/D  ■  1.5,  "  0.9.  (Adapted  from  Ref.  2). 
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Realizing  this,  several  researchers  have  set  about  to  capture  the 
Important  nonlinear  aspects  of  the  high-speed  rotor  problem  (Refs.  1,  11-13) 
Two  different  but  fundamentally  similar  approaches  have  been  taken: 

(1)  nonlinear  "quadrupole"  terms  have  been  re- Incorporated  Into  the 
"acoustic  analogy"  formulation  as  additional  source  terms;  and  (2)  a  non¬ 
linear  potential  equation  has  been  formulated  for  the  acoustic  radiation 
of  rotors. 

In  this  paper,  we  compare  and  highlight  the  differences  and  similari¬ 
ties  of  both  methods.  The  potential  formulation  will  first  be  used  to  help 
explain  some  of  the  more  Interesting  experimental  features  of  the  nonlinear 
flow  field  surrounding  the  blade.  Then  the  acoustic  analogy  approach  will 
be  discussed  and  quadrupole  sources  will  be  used  to  help  predict  the  result¬ 
ing  noise  field.  Finally,  the  results  will  be  critically  compared  with 
specially  run  tests  to  highlight  the  promise  or  fallings  of  the  acoustic 
analogy  approach  and  to  explain,  in  simple  physical  terms,  how  sound  Is 
radiated  by  a  rotor  with  transonic  tip  velocities. 


2.  THE  POTENTIAL  EQUATKW  AND  ITS  IMPLICATIONS 

The  acoustic  radiation  of  an  isolated  rotor  hovering  at  transonic 
tip  Mach  numbers  Is  considered  in  this  paper.  The  more  difficult  problem 
of  high-speed  forward  flight  Is  deferred  until  the  steady  transonic  and 
acoustic  features  of  the  hovering  problem  have  been  Identified.  Assuming 
constant  specific  heats  and  weak  shocks  (l.e.,  negligible  entropy  Increases) 
the  classical  potential  equation  Is: 
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(1) 


4.  2  grad^  •  grad 

+  ^  grad^  •  grad[(grad^)^]  ■  0 

where  ^  represents  the  velocity  potential  and  a  is  the  local  speed  of 
sound.  This  nonlinear,  time-dependent  equation  is  written  in  a  space-fixed 
coordinate  system;  however,  the  aerodynamics  of  a  hovering  rotor  are  basi¬ 
cally  steady  tdien  viewed  from  a  blade-fixed  reference  frame.  Therefore, 
following  the  work  of  Isom  (Ref.  11)  the  governing  potential  equation  can 

cylindrical  coordinates  and  expanded  to 


♦ee  ' 

■  {“5  +  -  «»*e)  (♦„  *T*  ♦z.) 


ag  -  undisturbed  speed  of  sound 
r  >  radial  distance  from  the  axis  of  the  cylindrical  coordinate  system 
y  -  ratio  of  specific  heats 

This  nonlinear  but  steady  second-order  partial  differential  equation  gov¬ 
erns  the  acoustic  as  well  as  the  transonic  behavior  of  the  hovering  rotor. 
At  the  present  time,  no  closed-form  solutions  to  this  equation  exist.  A 
procedure  adopted  by  some  researchers  is  to  numerically  solve  limited 
regions  of  the  aerodynamic  flow  field  (Refs.  14,  15).  Others  (Ref.  11) 
have  chosen  to  solve  the  nonlinear  acoustic  far-fleld,  using  weak-shock 
theory  (Ref.  16).  As  we  shall  see,  neither  is  a  completely  satisfactory 
solution,  for  the  nonlinear  aerodynamic  and  acoustic  fields  are  Interwoven. 

The  choice  of  a  cylindrical  coordinate  system  whose  axis  is  aligned 
with  the  rotor  is  sketched  in  Figure  2.  An  observer  riding  in  this  coordi¬ 
nate  system  sees  a  free-stream  velocity  that  increases  linearly  from  zero 
at  the  origin  to  ur  at  r.  As  indicated,  this  Increasing  free-stream 
velocity  continues  out  past  the  tip  of  the  rotor;  it  will  be  shown  to  be 
Important  to  many  of  the  arguments  to  come. 

Before  attempting  to  solve  Equation  (2) ,  it  is  instructive  to  follow 
the  approach  of  References  11  and  12  and  explore  the  behavior  of  the  gov¬ 
erning  equation.  It  is  known  from  the  theory  of  partial  differential 
equations  that  the  coefficient  of  4^^  governs  the  general  character  of 
the  potential  equation.  For  when 

A  H  0)2  -  ^  -  (Y  +  1)  ^  40  (3) 

<  0  :  elliptic  behavior 
>  0  :  hyperbolic  behavior 


be  transformed  to  blade-fixed 
second-order,  yielding 

|o.2  -  -  (Y  +  1)  ^  *e| 

where 

0)  "  angular  rotation  rate 


However,  A  takes  a  more  recognizable  form  after  some  further  manipulation: 


A 


(Y  -  1)(D 
.2 


Using  the  energy  equation,  this  becomes 
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r2 


(4) 


(5) 


Defining  U  ■  ur,  ■  -ur,  we  obtain 

09 
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(U»  +  u)^  u2) 


(6) 


The  last  term  In  the  brackets  of  Equation  (6)  can  be  neglected.  It  Is 
second  order  to  A  and  a  third-order  correction  to  the  second-order  poten¬ 
tial  equation.  If  we  now  define 


U„  +  u 

M,  = - -  local  Hach  number, 

1  a 

the  coefficient  of  (Eq.  (6))  becomes 

A  -  (1  -  mJ)  (7) 

He  have  shown  that  Che  general  behavior  of  the  second-order  transonic 
potential  equation  is  governed  by  the  local  Hach  nwtber  of  the  flow.  If 
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<  1.0,  Chen  A  <  0  and  Che  governing  equaclon  Is  elllpClc.  In  chls 
case,  no  wave-llke  scruccure  Is  possible.  However,  If  >  1.0,  Chen 
A  >  0  and  Che  governing  nonlinear  parclal  differencial  equaclon  Is  hyper¬ 
bolic.  Then,  characcerlsClcs  are  formed  along  which  dlsCurbances  can  propa— 
gaCe  In  a  wave-llke  manner.  IC  Is  also  ImporCanc  Co  realise  ChaC  Che  local 
Mach  nuiiri>er  Is  dependenc  on  Che  free-sCream  velocicy  ur,  Che  local 
speed  of  sound  a,  and  Che  local  percurbaclon  velocity,  u  >  -^g/r.  Each 
conCrlbuClon  Is  separaCely  idenclflable  in  Equaclon  (6). 

These  Ideas  are  quice  useful  when  one  is  accempclng  Co  InCerpreC  the 
experimental  results  of  Reference  2.  In  chat  sequence  of  experiments,  a 
1/7-scale  twisted  and  untwisted  UH-IH  model  rotor  (NACA  0012  airfoil)  was 
run  In  an  anecholc  hover  chsmber  at  near  zero  thrust  at  high  subsonic  tip 
Mach  numbers.  In-plane  acoustic  signatures  were  measured  and  stationary 
hot-wire  measurements  of  the  local  perturbation  field  (u)  surrounding  Che 
blade  were  made.  A  phenomenon  labeled  "delocalization"  was  observed  In 
this  experiment  as  the  hover  free-stream  tip  Mach  number  was  Increased 
above  0.88.  Shock  waves  on  the  rotor  suddenly  "delocalized"  and  radiated 
to  the  acoustic  far-fleld.  The  measured  far-fleld  acoustic  signature 
transitioned  from  a  near  symmetrical  triangular  pulse  shape  to  a  distinct 
sawtooth  character  —  confirming  the  existence  of  a  radiating  shock  wave 
(Fig.  3). 

The  connection  between  the  phenomenon  of  delocalization  and  the 
change  In  equation  type  has  been  suggested  In  References  11  and  13  and 
reported  In  Reference  12,  based  on  a  numerically  calculated  local  tran¬ 
sonic  flow  field.  In  the  following  paragraphs,  the  relationships  are 
experimentally  confirmed  and  are  shown  to  depend  on  the  local  Mach  nuiid>er 
of  the  flow.  Three  distinct  cases  are  considered:  ■  0.85, 

M^j^P  ■  0.88,  and  ■  0.90.  Some  degree  of  graphical  freedom  is  taken 

In  the  Interest  of  presenting  a  clear  sketch  of  the  basic  relationships 
Involved.  The  data  are  the  same  as  those  reported  In  Reference  17,  but  are 
expanded  upon  here.  In  the  figures  that  follow,  the  top  views  are  sketches 
of  events  pieced  together  with  limited  experimental  data,  and  the  aft  views 
are,  for  the  most  part,  Interpolations  of  experimental  data. 


Mt  -  0.M  M-f  -  0.89  My  -  0.90 

(s)  (b)  (e) 


flginr*  3>  Waveform  transition;  the  development  of  a  radiating  discontinuity 
la-plans,  r/O  -  1.5,  M^.^  -  0.9. 


Figure  4  depicts  the  top  and  aft  views  of  a  rotor  operating  at  a 
free-stream  tip  Mach  number  of  0.85.  A  locally  supersonic  region 

exists  near  the  tip  of  the  rotor.  For  this  region  M£  ■  (ur  +  u)/a  >  1.0, 
even  though  wr/sQ  all  along  the  blade  span  Is  less  than  0.85.  The  hyper¬ 
bolic  nature  of  this  pocket  of  supersonic  flow  is  a  result  of  local  aero- 
d3mamlc  nonllnearltles  (l.e.,  changes  in  the  local  speed  of  sound  a,  and 
local  perturbation  velocity,  u).  Surrounding  this  locally  supersonic  flow 
is  a  subsonic  flow  region,  Mj^  <  1.0,  in  which  the  governing  potential 
equation  is  elliptic.  No  wave-like  behavior  is  possible  through  this  com¬ 
pressible  elliptic  region.  However,  as  r  increases  beyond  the  tip  of  the 
blade,  again  becomes  greater  than  1  because  of  the  linearly  increasing 
free-stream  velocity  field  of  the  blade-fixed  cylindrical  coordinate  system 
For  this  region  u  «  0,  a  «  sq,  so 

M«  ~  >  1.0 
*  ao 

The  surface  where  this  first  happens  has  been  called  the  sonic  cylinder 
(Refs.  11,  13).  At  radii  larger  than  the  sonic  cylinder,  the  equation 
again  becomes  hyperbolic,  and  wave-like  propagation  is  certain.  The  acous¬ 
tic  implications  of  this  0.85  case  begin  in  the  hyperbolic  pocket  of  flow 
near  the  blade  tip.  Wave-like  disturbances  in  this  region  terminate  on 
the  boundary  of  an  elliptic  region  where  they  no  longer  propagate  in  char¬ 
acteristic  directions.  The  wave-like  character  of  the  inner  podket  is  thus 
broadened  as  information  passes  through  the  elliptic  region  to  the  sonic 
cylinder.  These  broadened  disturbances  are  then  propagated  in  a  wave-like 
manner  throughout  Che  outer  hyperbolic  region.  The  result  is  a  smoothly 
varying  near  symmetrical  acoustic  signature  in  the  far-fleld. 
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Figure  4.  Top  and  aft  views  of  shock  boundaries 


of  a  rotor,  M. 


0.85. 


The  competing  phenomena  become  even  more  interesting  when  the  undis¬ 
turbed  free-stream  tip  Mach  number  of  the  rotor  (M^ip)  is  increased  to  0.88 
(Fig.  5).  The  inner  supersonic  (or  hyperbolic)  region  grows  and  extends 
off  Che  tip  of  the  rotor  again  being  driven  by  local  aerodynamic  non- 
linearities.  At  the  same  time,  Che  higher  free-stream  tip  Mach  number  of 


LINEAR  SONIC 


Figure  S.  Top  and  aft  views  of  shock  boundaries  of 
a  rotor,  ■  0.88. 


the  rotor  decreases  the  radius  of  the  sonic  cylinder,  thus  moving  the  outer 
hyperbolic  region  toward  the  rotor  tip.  In  addition,  the  proximity  of  the 
linear  sonic  cylinder  to  the  blade  tip  introduces  aerodynamic  nonlinearities. 
These  tend  to  warp  the  sonic  cylinder  Inward,  bringing  the  two  hjrperbolic 
regions  even  closer  together.  However,  the  hyperbolic  regions  do  not  over¬ 
lap,  thus  Insuring  that  locally  generated  waves  In  the  Inner  region  do  not 
propagate  along  characteristics  to  tue  hyperbolic  far-field.  '  Instead,  the 
wave-llke  disturbances  are  forced  to  pass  through  a  small  elliptic  region 
where  they  are  broadened  before  entering  the  outer  hyperbolic  region  for 
propagation  to  the  far-fleld.  The  resulting  acoustic  signature  becomes 
more  sawtoothed  In  character,  but  does  not  contain  radiating  shocks. 

The  last  and  most  Interesting  condition,  in  which  the  local  free- 
stream  tip  Mach  number  Is  Increased  to  0.9,  Is  sketched  In  Figure  6.  The 

LINEAR  SONIC 
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localized  inner  hyperbolic  and  outer  hyperbolic  regions  connect  off  the 
blade  tip,  forming  one  continuous  supersonic  region  (M^  >  1.0).  In  this 
case,  shock  waves  that  are  generated  on  the  surface  of  the  rotor  now  propa¬ 
gate  uninterruptedly  In  a  radially  outward  characteristic  direction  to  the 
acoustic  far-fleld.  The  resulting  delocalization  phenomenon  Is  quite 
striking,  for  the  character  and  the  Intensity  of  the  acoustic  signature 
change  dramatically.  At  all  three  of  these  conditions,  measured  values  of 
local  Mach  number  support  and  explain  the  phenomenon  of  "transonic 
delocalization. " 


3.  THE  QUADRUPOLE  APPROXIMATION  AND  THE  ACOUSTIC  FIELD 

Although  the  phenomenon  of  delocalization  has  been  explained  by 
simply  looking  at  the  coefficient  of  4)00  in  Equation  (2) ,  predicting  the 
radiating  acoustic  field  Is  another  matter.  The  explanations  presented  are 
themselves  functions  of  either  measured  or  calculated  flow  properties.  In 
essence,  a  near-fleld  description  of  the  aerodynamic  flow  field  Is  required 
before  the  events  In  the  acoustic  far-fleld  can  be  explained.  Even  then, 
acoustic  nonllnearltles  in  Equation  (2)  may  alter  the  waveform  of  the  prop¬ 
agating  wave  (Refs.  11,  13).  Precise  calculations  of  the  radiating  sound 
field  are  dependent  on  the  full  solution  of  the  nonlinear  potential  equa¬ 
tion  (Eq.  (2)). 

On  the  other  hand,  the  successful  explanation  of  the  delocalization 
phenomenon  suggests  that  local  aerodynamic  nonllnearltles  strongly  Influ¬ 
ence  the  acoustic  radiation  problem.  Therefore,  a  logical  step  In  the 
calculation  of  the  acoustic  field  Is  the  Incorporation  of  the  near-fleld 
aerodynamic  nonllnearltles  In  the  acoustic  radiation  equation. 

At  least  two  alternative  ways  of  Implementing  this  Idea  have  been 
presented  In  the  literature.  One  recently  proposed  method  (Ref.  12)  maps 
the  nonlinear  near-fleld  to  a  nonrotating  control  surface  where  Klrchoff's 
theorem  Is  applied.  As  reported,  the  control  surface  must  be  chosen  to  be 
large  enough  to  capture  the  nonlinear  aerodynamic  behavior  of  the  problem, 
but  not  to  be  so  large  as  to  make  numerical  computation  Impractical.  Ini¬ 
tial  calculations,  using  this  procedure  coupled  with  an  existing  near-fleld 
numerical  code,  have  shown  Improvement  In  peak  amplitude  levels,  but  not 
much  Improvement  In  waveform  characteristics  above  the  delocalization  Mach 
number.  As  discussed  In  Reference  12,  this  Is  most  likely  due  to  the 
numerical  Insensitivity  of  the  transonic  code  at  the  boundary  of  the  non¬ 
rotating  control  surface. 

The  second  method  employs  the  well-known  "acoustic  analogy"  proce¬ 
dures  to  evaluate  the  volume  distributions  of  local  aerodynamic  nonllnearl¬ 
tles  or  quadrupoles  (Refs.  1,  5,  17).  This  approach  Is  utilized  in  this 
paper  together  with  an  accurate  estimation  of  the  local  three-dimensional 
nonlinear  aerodynamic  flow  field.  It  was  first  attempted  In  Reference  1 
for  high-speed  helicopter  Impulsive  noise  by  approximating  the  quadrupole 
term,  using  two-dimensional  Frandtl-Glauert  aerodynamic  theory.  Although 
small  Improvements  In  amplitude  were  noted,  there  remained  major  discrep¬ 
ancies  between  theory  and  experiment.  The  authors  of  Reference  5  performed 
a  similar  calculation  for  a  transonic  propeller  using  two-dimensional 
transonic  codes  for  Input;  they  reported  better  agreement  In  amplitude  but 
poor  agreement  In  pulse  shape.  Three-dimensional  transonic  aerodynamics 
vere  Included  In  Reference  17  for  the  hovering  helicopter  rotor,  and  the 
results  were  more  Impressive.  However,  the  limitation  of  Integration 
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region  for  the  quadrupole  term  and  the  omission  of  shock  terms  caused  poor 
agreement  between  pulse  shape  and  theory  near  or  above  the  delocalization 
Mach  number.  In  this  paper,  these  latter  two  restrictions  are  removed 
from  the  quadrupole  approximation  In  the  transonic  acoustic  field. 


The  analysis  begins  with  the  well-known  Ffowcs  Williams  and  Hawkings 
formulation  (Ref.  7),  In  which  the  sound  radiated  by  surfaces  In  motion  Is 
expressed  as  the  Integral  equation: 


Airagp 


32  fr 

3x .  3x ,  RI  1  - 
1  1  L  ' 
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(8) 


where 


"il  - 


and  standard  notation  Is  assumed. 


The  first  term  In  Equation  (8)  Is  the  linear  thickness  contribution 
to  the  radiated  noise.  This  surface  Integral  has  been  evaluated  by  many 
authors,  but  their  efforts  yielded  the  generally  poor  agreement  with  exper¬ 
iment  discussed  previously.  The  second  term,  also  a  surface  Integration 
over  the  blade,  can  be  classified  as  either  a  linear  or  a  nonlinear  term, 
depending  on  how  the  surface  pressure  Is  approximated.  However,  for  most 
of  the  in-plane  computations  reported  here,  its  contribution  to  the  radi¬ 
ated  noise  field  Is  small  and  thus  has  been  neglected.  A  strict  evaluation 
of  the  quadrupole  (third)  term  in  Equation  (8)  requires  a  volume  integra¬ 
tion  over  all  space.  However,  the  Interest  In  this  paper  Is  In  capturing 
those  nonlinear  aerodynamic  terms  that  may  govern  the  acoustic  radiation 
problem.  Therefore  the  evaluation  of  the  quadrupole  term  Is  confined  to  a 
volume  Integration  within  the  aerodynamic  nonlinear  flow  field  of  the  blade. 


The  quadrupole  term  assume:;  a  simpler  form  If  we  restrict  our 
attention  to  the  acoustic  far-field.  Then  the  spatial  differentiations 
can  be  easily  converted  to  time  differentiations.  The  last  term  of 
Equation  (8)  becomes 


3^ 

3x^3x 


j 


(9) 


where  Tjq^  ■  '^lj^*^j>  8  is  the  unit  vector  from  the  source  at  the 

retarded  time  to  an  observer  In  the  acoustic  far-fleld.  It  Is  known  from 
transonic  computations  and  experimentation  that  the  primary  quadrupole 
regions  are  confined  within  a  few  chord  lengths  normal  to  the  rotor  plane. 

For  In-plane  far-fleld  radiation,  the  unit  vector  R  Is  nearly  in 
the  blade  rotational  plane  and  Is  nearly  parallel  to  the  blade  chordwlse 
direction  when  the  acoustic  pressure  reaches  Its  peak  level.  If  Isentroplc 
flow  Is  assumed  and  the  perturbation  velocities  are  measured  In  the  coordi¬ 
nate  system  given  In  Figure  2,  T{^  becomes 


1 
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Pnfv?  cos^6  +  2v  V.  cos9  slnS  +  sin^S 
0\  e  re  r  , 


(10) 


where  the  z  component  of  the  perturbation  velocity  does  not  appear 
because  of  the  choice  of  an  in-plane  far-field  microphone  position.  For 
simplicity  in  the  resulting  calculations,  it  has  also  been  assumed  that 
sin  0  s*  0  and  that  u  =«  V0  near  the  integration  region  of  interest. 
This  is  true  as  long  as  the  quadrupole  field  is  in  fact  localized  to  a 
region  near  the  rotor  tip.  Then,  Equation  (10)  becomes 

^RR  -  Po(^) 


where  u  represents  the  perturbation  velocity  along  the  blade  chord  and 
wr  is  the  free-stream  velocity  of  the  point  in  the  flow  field  being  eval¬ 
uated.  The  two  terms  represented  in  Equation  (10)  arise  from  similar 
properties  of  the  flow  already  discussed  in  the  potential  formulation 
(Eq.  (4)).  Changes  in  the  local  speed  of  sound  and  local  streamwlse  per¬ 
turbation  nonlinearities  are  included  although  the  equation  forms  do  not 
permit  a  one-to-one  correspondence  of  terms. 

Equations  (8),  (9),  and  (11)  describe  the  nonlinear  far-field  acous¬ 
tic  radiation  of  the  transonic  hovering  rotor.  For  subsonic  tip  Mach 
numbers,  numerical  evaluation  of  the  surface  integrals  presents  no  real 
problems.  However,  the  volume  integration  of  quadrupoles  is  not  as 
straightforward. 

As  discussed  previously,  the  accuracy  of  this  nonlinear  formulation 
is  totally  dependent  on  the  detailed  knowledge  of  the  flow  field  surround¬ 
ing  the  blade.  In  this  paper,  the  flow  fields  of  interest  were  computed 
by  the  three-dimensional  transonic  numerical  code  described  in  Reference  15. 
This  code  solves  the  near-field  transonic  small-disturbance  potential  equa¬ 
tion  in  a  blade-fixed  reference  frame.  The  solution  uses  a  conservative, 
mixed-difference  relaxation  scheme.  Because  this  code  was  developed  to 
predict  blade-surface  aerodynamics,  some  simplifying  assumptions  were  made 
to  Equation  (2)  —  namely,  the  nonlinearities  known  to  be  small  near  the 
blade  surface  were  neglected.  In  particular,  the  terms  2u)()>^(|>^0  and 
2(«><^z<fz0  0*1  left-hand  side  and  the  factor  (y  ~  l)u)(()0  on  the  right-hand 

side  of  Equation  (2)  were  not  programmed. 

The  nature  of  the  calculated  chordwise  pressure  distributions  as  a 
function  of  blade  radius  is  shown  in  Figure  7  for  a  tip  Mach  number  of  0.9 
at  zero  angle  of  attack.  On  the  blade  surface,  a  shock  wave  forms  at 
about  90%  of  the  span  and  persists  over  the  outer  portions  of  the  blade 
and  beyond.  These  calculated  pressure  distributions  are  consistent  with 
the  experimentally  sketched  flow  regions  of  Figure  6.  The  maximum  value 
of  shock  strength  is  calculated  to  occur  at  95%  radius.  It  is  also  impor¬ 
tant  to  note  that  the  shocks  that  are  calculated  with  this  code  do  not  look 
like  measured  shock  waves.  The  numerical  solution  has  built  into  it  a 
numerical  viscosity  which  tends  to  smooth  discontinuities  over  several  mesh 
points.  This  tends  to  smooth  the  final  acoustic  waveforms  to  some  degree. 


Cp  ON  THE  BLADE  SURFACE 


- Cp  AT  1.18  CHORD  UP  VERTICALLY 
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1.068 


V3/Ro  =  1.038 


TIP  MACH  NUMBER  =  0.9 
ASPECT  RATIO  =  13.71 
ZERO  INCIDENCE 

Figure  7.  Cp  distribution  on  the  blade  surface  and  at 
1.18  chords  up,  aspect  ratio  -  13.7,  NACA  0012,  near  zero 
angle  of  attack,  M  -  0.9. 


The  vertical  extent  of  the  local  flow  field  is  shown  in  Figure  8  for 
the  95%  span  location.  The  rapid  decay  in  disturbance  velocity  u  has 
been  confirmed  by  experiment.  However,  as  indicated  in  Reference  17,  the 
decay  with  vertical  distance  from  the  airfoil  is  much  faster  than  with 
two-dimensional  transonic  theory,  but  much  slower  than  with  linear  calcu¬ 
lations.  The  faster  decay  of  three-dimensional  transonlcs  is  caused  by 
tip  relief.  This  figure  illustrates  that  the  blade  tip  region  is  governed 
by  local  nonlinear  three-dimensional  transonic  efforts  and  must  be  calcu¬ 
lated  accordingly. 


MAX.CHORDWISE  PERTURBATION 
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Figure  8.  Vertical  extent  of  the  local  flow 


field. 
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The  final  problem  in  the  evaluation  of  the  quadrupole  Integral  is 
the  development  of  a  calculation  procedure  for  Equations  (8) ,  (9) ,  and  (11) 
that  is  valid  when  Mr  =  1.0.  This  problem  occurs  when  the  volume  integra¬ 
tion  is  extended  up  to  and  beyond  the  linearized  sonic  cylinder.  The  inte¬ 
grand  in  Equation  (9)  then  contains  the  product  of  two  terms  which  compete 
to  decide  the  eventual  magnitude  of  the  quadrupole  radiation.  The  first 
is  the  decaying  source  field  represented  by  the  TpR  term  in  Equation  (9) . 
This  is  multiplied  by  the  1/ | 1  -  Mr|  term  which  goes  to  “  as  Mr 
approaches  1.  Fortunately,  the  singularity  is  Integrable,  but  it  must  be 
handled  quite  carefully.  In  this  paper  the  acoustic  planform  technique 
was  chosen  to  perform  the  numerical  integration  near  Mr  ->  1.0.  A  complete 
discussion  of  the  procedures  and  pertinent  references  is  given  in  the 
Appendix. 

Equation  (11)  governs  the  magnitude  of  the  quadrupole  source  strength, 
Trr.  For  an  observer  in  the  acoustic  far-field,  and  for  an  integration 
region  confined  near  the  blade  tip,  all  effective  source  points  normal  to 
the  rotor  plane  will  arrive  at  the  same  retarded  time.  Therefore  it  is 
possible  to  Integrate  u^  in  Che  z  direction  and  effectively  collapse 
the  quadrupole  volume  Integration  to  a  surface  Integral  in  the  rotor  plane 
(Ref.  17).  The  resulting  integration  for  the  model  rotor  operating  at  near 
zero  thrust  at  a  hover  tip  Mach  number  of  0.9  is  shown  in  Figure  9  for  sev¬ 
eral  radial  stations.  The  integral  dz  is  the  largest  at  about 

r/ro  =  0.95  and  decreases  both  inboard  and  outboard  of  this  point.  The 
monotonlc  decrease  in  amplitude  out  past  the  tip  to  beyond  the  linear  sonic 
circle  suggests  that  quadrupole  sources  die  quickly.  However,  the  ultimate 
strength  of  the  quadrupole  term  is  dependent  on  the  integration  over  an 
effective  planform  of  the  product  of  J\i^  dz  and  l/|l  -  Mr| .  Because 
l/|l  -  MrI  is  singular  when  Mr  >  1.0,  the  net  quadrupole  Integration  is 
weighted  heavily  near  the  sonic  circle.  The  resulting  contributions  to  the 
acoustic  field  are  shown  in  Figure  10  for  the  same  hovering  rotor  at  near 
zero  thrust.  As  shown,  the  dominant  integration  region  for  the  quadrupole 
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Figure  9.  Quadrapole  source  strength  as  a  function  of  r/r^ 
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term  Is  between  the  rotor  blade  tip  and  the  linear  sonic  circle.  The  figure 
also  Illustrates  how  the  shape  of  the  quadrupole  term  changes  near  a  hover 
tip  Mach  number  of  0.9  for  the  UH-IH  model  rotor.  At  these  high  tip  speeds, 
the  shape  of  the  quadrupole  term  Is  strongly  Influenced  by  the  shape  of 
Ju^  dz  near  the  sonic  circle.  The  resulting  waveform  becomes  sawtoothed 
In  character  because  the  governing  equations  are  nearly  hyperbolic  In  char¬ 


acter  In  this  region.  The  structure  of  the  Ju^  dz  field  Is  nearly  pre¬ 


served,  yielding  a  shock-llke  quadrupole  waveform  In  the  acoustic  far-fleld. 


This  does  not  happen  at  lower  hover  tip  Mach  numbers.  The  amplitude 
of  the  Ju^  dz  field  dies  away  rapidly  causing  the  quadrupole  Integral  to 
be  dominated  by  the  linear  elliptic  region.  Here,  no  wave-llke  structure 
Is  possible  and  the  structure  of  Ju^  dz  Is  smoothed  In  the  amplification 
process,  producing  nearly  synmetrlcal  pressure  time  histories. 
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Figure  10.  Radial  cootrl 
butlon  of  quadrapole 
radiation. 


COMPARISON  WITH  EXPERIMENT 


An  evaluation  of  the  overall  prediction  accuracy  is  presented  below 
by  a  comparison  with  data  presented  in  Reference  2.  In  that  test,  a 
1/7  scale  model  of  a  UH-IH  helicopter  (NACA  0012  airfoil  with  10.9"  twist 
and  with  no  twist)  was  run  at  nearly  zero  thrust  at  transonic  tip  Mach 
numbers  from  0.8  to  1.0  in  hover  in  an  anecholc  hover  chamber.  Micro¬ 
phones  located  in  the  plane  of  the  rotor  recorded  the  radiated  noise. 

Figure  11  presents  the  monopole  and  quadrupole  contributions  to  the 
radiated  noise  at  M^j^p  =  0.8.  As  discussed  previously,  at  these  low  tip 
Mach  nuDd>ers  the  monopole  contribution  to  the  radiated  In-plane  far-fleld 
noise  underpredicts  the  magnitude  of  the  measured  data  (Ref.  2),  but 
accurately  predicts  the  pulse  shape  (Fig.  11a) .  The  quadrupole  contribu¬ 
tion  is  of  an  almost  similar  shape  (Fig.  11b)  and  helps  Improve  the  corre¬ 
lation  with  experiment  (Figs.  11c,  lid)  to  nearly  acceptable  levels. 

As  the  hover  tip  Mach  number  is  increased  to  0.88  (Fig.  12)  the 
monopole  and  quadrupole  terms  (Figs.  12a,  12b)  become  equally  large  in 
magnitude.  More  Importantly,  however,  the  shape  of  the  quadrupole  term 
begins  to  change  from  a  near-symmetrical  to  an  asymmetrical  character. 

This  is  noted  on  the  pressure  recovery  side  of  the  quadrupole  calculation 


and  Is  a  direct  result  of  Including  the  local  shock  structure  of  the  flow 
field  In  the  numerical  calculation.  When  the  monopole  and  quadrupole 
contributions  are  added,  good  correlation  In  amplitude  and  pulse  shape  Is 
observed  (Figs.  12c,  12d).  The  overall  shape  of  theory  and  experiment  are 
still  basically  symmetrical  In  character,  but  the  local  shock  structure  of 
the  transonic  flow  field  Is  acting  to  destroy  this  symmetry. 

At  the  slightly  higher  hover  tip  Mach  number  of  0.90,  localized 
transonic  effects  cause  large  changes  to  the  radiated  noise  field  (Fig.  13) 
Although  the  linear  term  (monopole.  Fig.  13a)  remains  quite  symmetrical  In 
shape  and  substantially  underpredicts  the  measured  data,  the  nonlinear 
(quadrupole.  Fig.  13b)  term  changes  shape  dramatically  and  Increases  In 
amplitude.  This  is  a  reflection  of  the  fact  that  local  shocks  are  propa¬ 
gating  to  the  far  acoustic  field  ("delocalization”  In  Ref.  17).  When 
the  monopole  and  quadrupole  terms  are  compared  with  experimental  data 
(Figs .  13c ,  13d) ,  good  agreement  In  pulse  shape  Is  observed .  To  the 
authors'  knowledge,  this  Is  the  first  time  calculated  pulse  shapes  have 
compared  favorably  with  measured  data  for  a  rotor  operating  near  the 
"delocalization"  Mach  number. 

The  relative  accuracy  with  which  the  peak  negative  amplitude  of  the 
high-speed  Impulsive  noise  phenomenon  can  be  predicted  Is  Illustrated  In 
Figure  14.  Much  better  agreement  between  theory  and  experiment  Is  demon¬ 
strate!  for  Mach  numbers  of  up  to  0.9.  At  M^ip  -  0.9,  amplitudes  are 
overpredicted.  This  result  is  related  to  the  sensitivity  of  this  approach 
to  perturbation  velocities  near  the  linear  sonic  circle  and  to  the  accuracy 
with  which  existing  transonic  near-field  codes  can  predict  these  velocities 
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Figure  14.  Comparison  of  theory  and  experiment  In 
hover,  In-plane,  r/D  ■>  1.5. 


5.  SOME  FDRTHER  OBSERVATIONS 


Some  insight  into  why  overprediction  of  the  acoustic  pulse  near  the 
"delocalization”  Mach  number  is  likely  can  be  gained  from  Figure  15.  In 
this  figure,  peak  values  of  measured  horizontal  perturbation  velocities 
are  compared  with  the  local  transonic  potential  code  used  in  this  paper 
for  Mfip  ■  0.9.  Inward  from  the  tip  of  the  rotor  blade  comparison  with 
experiment  is  quite  good;  however,  near  the  linear  sonic  cylinder  where 
the  quadrupole  contributions  to  the  radiated  noise  are  large,  theory  and 
experiment  diverge.  It  is  quite  likely  that  the  overestimation  of  the 
transonic  perturbation  velocity  in  this  region  causes  the  quadrupole  con¬ 
tribution  (which  squares  the  error)  to  become  too  large.  A  transonic  code 
that  incorporates  all  of  the  near-field  nonlinearities  should  help  the 
correlation  of  this  acoustic  theory  with  experiment  at  the  higher  transonic 
tip  speeds. 


A  plot  of  the  torque  coefficient  versus  hover  tip  Mach  number  is 
given  in  Figure  16  for  the  untwisted  model  rotor  (Ref.  2)  at  near  sero 
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thrust.  The  resulting  curve  Increases  sharply,  defining  a  point  at  a  hover 
tip  Mach  number  of  about  0.85  that  might  be  considered  to  be  the  onset  of 
"drag  divergence."  This  Is  in  contrast  to  the  phenomenon  of  "delocaliza¬ 
tion,"  which  occurs  at  >  0.9  for  this  model  rotor.  It  Is  probable 

that  the  two  events  are  related,  but  the  experimental  and  theoretical  evi¬ 
dence  to  date  does  not  explicitly  connect  drag  rise  with  radiated  noise. 


6.  CONCLUSIONS 


Inclusion  of  the  quadrupole  terms  In  the  numerical  evaluation  of  the 
Ffowcs  Williams  and  Hawkings  Integral  equation  improves  the  comparison  with 
measured  data  for  high-speed  helicopter  rotors  operating  at  tip  Mach  num¬ 
bers  from  0.8  to  0.9.  The  predicted  waveform  closely  matches  the  measured 
data.  At  low-tip  Mach  numbers  (M^ip  ~  0.8),  nearly  symmetrical  pulse  shapes 
are  predicted  and  measured  while  at  the  delocalization  Mach  number 
(Mtip  ~  0.9)  sawtoothed  shock-like  waveforms  are  predicted  and  measured. 
However,  at  Mach  numbers  greater  than  0.9,  theory  tends  to  overpredict  mea¬ 
sured  peak  negative  amplitudes. 

The  changing  character  of  the  acoustic  waveform  depends  critically 
on  the  local  three-dimensional  transonic  flow  field  surrounding  the  blade 
tip.  If  the  local  supersonic  region  surrounding  the  blade  Is  bounded  by 
a  subsonic  cylindrical  region,  the  waveform  remains  near  symmetric  In 
character.  However,  at  higher  Mach  numbers  (Mt-j^p  ■  0.9),  when  the  local 
supersonic  region  connects  with  an  outer  hyperbolic  region,  the  waveform 
becomes  sawtoothed  and  local  shocks  propagate  (delocalize)  to  the  acoustic 
far-f leld . 

Predicted  acoustic  waveforms  are  dependent  upon  accurate  knowledge 
of  the  near-fleld  perturbation  velocity  field  near  the  blade  tip.  For  the 
present  paper,  this  near-field  Is  computed  using  a  small  disturbance  code 
which  Incorporates  near-fleld  aerodynamic  nonllnearltles  near  the  blade 
surface.  Measurements  of  peak  perturbation  velocities  indicate  that  this 
small  disturbance  transonic  theory  Is  in  good  agreement  In-board  of  the 
tip  region  but  becomes  less  accurate  (overpredicts)  off  the  blade  tip. 

This  Is  thought  to  be  one  reason  why  the  quadrupole  acoustic  formulation 
tends  to  overpredict  measured  values  at  Mach  nuiid>ers  greater  than  0.9.  It 
Is  doubtful  that  quantitative  prediction  of  far-fleld  transonic  rotor  noise 
at  Mach  numbers  greater  than  0.9  can  be  attempted  until  the  transonic  poten¬ 
tial  equation  Is  solved  over  regions  of  larger  radial  extent  than  the  equa¬ 
tion  has  dealt  with  to  date.  It  Is  clear  that  local  transonic  aerodynai^cs 
strongly  Influence  the  high-speed  rotor's  acoustic  far-fleld.  In  the  next 
few  years,  design  changes  based  on  this  new  understanding  should  help 
reduce  the  radiated  noise. 


APPENDIX 


The  Integral  equation  for  the  volume  displacement  effect  Is  shown 
as  follows  (Ref.  7): 


4ita|p'(8,t) 


it  Iff  |Vf  |6(g)dy  dx 


(Al) 


ifhere 


Q(y.T)  - 


T  ~  retarded  tine 
c  ■  observer  tine 
f  *  0:  body  surface 


(A2) 

R  -  |x  -  y| 

X  ■  observer  position  vector 
y  -  source  position  vector 


Since  the  Intervals  of  the  integrals  are  infinite,  the  order  of 
integration  can  be  interchanged.  Then  Equation  (Al)  becomes 

4ira2p*(5E,t)  -  6{f (y.x)}  |7f  (y.r)  |dy  (A3) 

where 


T  “ 
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If  the  7  coordinate  system  is  changed  to  the  1  coordinate  system,  such 
that  the  axis  is  perpendicular  to  the  surface  f(y,T)  >  0,  and  the 

CatCa  axes  are  on  that  same  surface  (see  sketch).  Equation  (A3)  becomes 


4ira2p'(x,t) 


-  —  jjj  (?.T)|  |Vf  (?,T)  I  J(y,?)d5  (A4) 


The  Jacobian  J(y,f)  between  coordinate  systems  y  and  ^  is  expressed 
as  follows 


irtiere 


This  expression  is  obtained  as  follows.  Define  a  function  7  such  that 


Then 


Also, 


|7F|  -  |7^f 


-  lylc* 

Therefore, 
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Now  Equation  (A4)  becomes  as  follows  with  the  Jacobian  J(y,C) 

JJJ"  qf(E.,)(d? 


(A5) 


Since  the  coordinate  system  Is  defined  such  that  Is  perpendicular  to 

the  surface  defined  by  f(f,T)  ■  0,  dC  becomes  as  follows: 

dC  -  dCi  d52  dCj 


“  df  dC2  dC3 

With  this.  Equation  (A5)  becomes 

4ira2p’(5E,t)  -  <*52  (A6) 

S: 

f^.T)  -  0 

This  Is  the  general  equation  for  the  acoustic  planform  method 
(derived  also  In  Ref.  8) ,  where  the  Integration  Interval  covers  only  the 
area  of  source  locations  In  retarded  time  for  which  signals  arrive  at  the 
observer  simultaneously.  This  is  equivalent  to  the  formulation  for  a 
collapsing  sphere  method  (see  Ref.  8).  The  area  S{f (0,S2»^3>''^)  " 
for  the  Integration  in  Eqxiation  (A6)  is  called  the  acoustic  planform. 
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The  acoustic  planform  changes  shape  depending  on  the  radiation 
(cooponent  of)  Mach  nuaiber,  t%.  For  very  low  )%,  the  acoustic  planform 
Is  almost  congruent  with  the  shape  of  the  blade  planform.  As  %  increases, 
the  acoustic  planform  distorts,  and  finally  breaks  into  multiple  regions  for 
%  greater  than  1,  as  shown  In  Figure  17.  Figure  17  shows  a  case  In  idiich 
a  blade  Is  rotating  In  a  circular  path  with  a  constant  supersonic  tip  Mach 
number  Mtip.  With  a  supersonic  tip  Mach  number,  the  radiation  component 
of  the  tip  Mach  number  varies  from  subsonic  to  supersonic,  and  then 
back  to  subsonic  speed  again  as  the  azimuthal  angle  Increases. 

The  multiple  regions  of  the  acoustic  planform  for  a  supersonic 
Is  explained  In  Figure  18. 


Mr 


Mr<1 


Mj 


Mr>1 


Mr<1 


Figure  17.  Physical  piste  and  acoustical  planform 
for  a  auparaonic  tip  Mach  nuMbar. 
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Figure  18.  Relative  position  between  wavefront 
and  blade. 


1.  At  t  0  a  sound  wave  is  generated;  because  the  speed  of  sound 
is  faster  than  the  blade  motion  Mji,  the  wavefront  moves  ahead  of  the  blade 
(see  Figs.  18a,  18b). 

2.  Mow,  the  blade  speeds  up  until  is  greater  than  1,  and  the 
blade  catches  up  with  the  wavefront  (Figs.  Idc-lde) . 
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3.  The  blade  again  slows  down  to  a  subsonic  speed,  and  the  wave- 
front  now  moves  ahead  of  the  blade  (Figs.  18f,  18g) . 

During  the  variation  of  the  blade  Mg  relative  to  the  propagation 
speed  of  the  wavefront,  the  wavefront  can  Interact  with  a  particular  point 
on  the  airfoil  one,  two,  or  three  times,  depending  on  the  value  of  the 
airfoil  Mr. 


These  multiple  values  of  source  positions  can  also  be  shown  from 
the  point  of  view  of  a  retarded  time  expression.  The  relationship  between 
the  observer  time  and  the  retarded  time  Is  given  as  follows: 

X  =  t--^ 

From  this  expression.  In  order  for  signals  from  two  separate  positions  on 
the  rotor  blade  to  reach  the  observer  simultaneously,  the  following  equa¬ 
tion  should  be  satisfied: 


where  and  R2  are  distances  from  two  different  sources  to  the  observer, 
and  i|)i  and  il(2  are  respective  azimuthal  angles  from  the  original  location; 
0)  Is  an  angular  velocity. 

For  a  given  angle  ipi,  the  position  of  the  second  source  \p2  may 
be  found  by  the  following  function: 

F(i|<)  =  (Ri  -  R2)  +  1/-2  -  'Pi  =  0  (A7) 

’ 


SOURCE  #1  POSITION 


SOURCE  #2 
POSITION  AT 
EMISSION 
TIME 


AT  EMISSION  TIME 


For  a  given  value  of  and  Rj ,  Figure  19  shows  the  graphical 
solution  for  Equation  (A7) .  That  is,  for  a  subsonic  Mach  number  Mr,  the 
function  F  Increases  monotonlcally  such  that  only  one  value  for  <|;2 
can  be  obtained  (see  Fig.  19a).  For  a  supersonic  case,  the  function  F 
has  one,  two,  or  three  roots,  depending  on  the  value  of  Mr;  that  Is, 
one  root  is  obtained  for  a  subsonic  Mr,  two  roots  for  Mr  ~  1,  and  three 
roots  are  obtained  for  Mr  >  1  (see  Fig.  19b) . 

A  typical  acoustic  planform  for  a  supersonic  tip  Mach  number  Is 
shown  in  Figure  20,  In  which  the  leading  and  trailing  edges  corresponding 
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Figure  20.  A  typical  acoustical  planfom  for  the 
caae  of  the  subsonic  tip  tfach  nuad>er. 


to  a  physical  blade  are  shown.  Two  observations  stand  out.  The  first  Is 
the  center  portion  of  the  planform,  where  the  trailing  edge  and  leading 
edge  are  reversed.  This  phenomenon  can  be  explained  from  Figure  18. 

When  a  wavefront  Is  ahead  of  an  airfoil  (see  Fig.  18c)  and  the  airfoil  Is 
accelerating,  the  wavefront  first  Interacts  with  the  leading  edge,  and 
then  sweeps  the  airfoil  to  the  trailing  edge  (see  Figs.  18d,  18e) . 

A  second  obseirvatlon  Is  that  of  the  two  local  depressions  A  and  B 
shown  In  Figure  20.  The  depth  of  these  depressions  depends  mainly  on  the 
magnitude  of  Mg.  Generally,  these  two  depths  have  different  values. 
Between  these  different  depths,  the  acoustic  planform  has  only  two  possible 
regions:  either  one  leading  edge  and  three  trailing  edges  or  three  trail¬ 
ing  edges  and  one  leading  edge,  depending  on  whether  Che  local  depres¬ 
sion  B  Is  deeper  or  shallower  than  A.  The  outermost  radial  portion  has 
three  regions,  and  the  Inner  radial  portion  has  one  region. 
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